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FOREWORD 


The  research  on  which  this  report  is  based  was  performed  by  the  author 
at  the  Cornell  Aeronautical  Laboratory,  Inc,  ,  Buffalo,  New  York,  under  Army 
Contract  DA  44- 177-TC-439,  Project  Number  9-38-01-000,  ST  902.  The 
Transportation  Corps,  U.  S.  Army  Transportation  Research  Command,  Fort 
Eustis,  Virginia,  is  the  monitoring  agency.  This  report  represents  part 
of  a  research  program,  which  is  devoted  to  the  investigation  of  several  specific 
problems  associated  with  STOL/VTOL  flight. 
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I. 


INTRODUCTION 


There  is  at  present  considerable  interest  in  the  aerodynamic  performance 

© 

of  the  vectored  slipstream  type  of  STOL  aircraft.  To  a  large  degree  the  per¬ 
formance  of  this  arrangement  depends  on  the  aerodynamic  interaction  of  two 
systems,  the  wing  and  the  propeller,  and  as  such  it  constitutes  a  formidable 
problem  in  flow  analysis.  In  the  interest  of  getting  wing  data,  one  simplifying 
approach  consists  of  neglecting  the  propeller,  per  se,  and  substituting  for 
its  system  a  jet  or  slipstream  of  indefinite  extent.  The  analysis  is  thereby 
reduced  to  one  of  an  isolated  wing  in  a  nonuniform  or  rotational  flow. 

It  is  usual  in  this  kind  of  analysis  to  effect  a  further  important  simplifi¬ 
cation  by  specifying  an  ideal  jet  and  postulating,  on  the  basis  of  small  disturbances, 
© 


that  the  cross-sectional  shape  will  remain  essentially  unchanged.  Under  such 
conditions  the  problem  is  related  in  some  respects  to  that  of  a  wing  in  an  open 
throat  wind  tunnel  -  -  a  scheme  that  has  been  examined  extensively  by  the 
method  of  images. 

The  latter  approach  has  been  applied  to  the  wing- jet  problem  by  a 
number  of  investigators,  whose  contributions  can  be  conveniently  grouped 

according  to  the  parameter,  jet  to  main- stream  velocity  ratio.  For  values 

,  ©  1 
of  this  parameter  close  to  unity,  the  problem  was  investigated  by  Koning  , 

2  3  4  5  6 

Franke  and  Weinig  ,  Ruden  ,  Bausch  ,  Kyewski  and  Vandrey  ,  Vandrey  , 

7  8  9 

etc.  For  large  values  we  have  the  developments  of  Glauert  ,  Stuper  ,  Pistolesi  , 

Baranoff^,  and  others.  The  more  recent  contributions,  namely  those  of 
11  12 

Rethorst  ,  and  Ribner  ,  seem  to  apply  for  all  values  of  the  jet-velocity  ratio 

13 

as  does  the  work  of  Graham  et  al 
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The  present  study  was  initiated  to  investigate  Jhe  effects  of  jet  curvature 

tjft 

on  the  premise  that  at  higher  angles  of  flow  deflection,  this  factor  should  have 

®  © 

an  important  influence  on  wing  lift.  The  jet  boundary  is  he^e  regarded  as  a 

vortical  layer,  as  is  also  the  wing  surface  and  its  wake.  The  jet  vorticity 
distribution  is  determined  subject  to  the  dynamic  equations  of  rotational  flow 
which  are  applied  in  a  thin  region  at  the  jet  boundary.  This  approach  leads  to 

(*£,  s  a, 

the  consideration  of  jet  trailing  vorticity  or  the-ea-called  -"secondary  vorticity" 

14  15 

which  occurs,  as  shown  by  Squire  and  Winter  and  also  by  Hawthorne  , 

Si  in  certain  curved  nonuniform  flows. 


AI-1190-A-6 


© 


II.  INTEGRAL  REPRESENTATION  OF  THE  GENERAL  VELOCITY  FIELD 

0 

To  show  clearly  the  present  approach,  it  will  be  convenient,  first,  to 
recall  certain  general  properties  of  vector  fields,  and  how  they  relate  to 
steady-flow  fields  of  an  incompressible,  inviscid  fluid  when  the  flow  field  includes 
free  vorticity.  We  begin  with  the  following  theorem  which  is  proved  in  Ref.  16 
Any  continuous  vector  field,  AJ  ,  defined  everywhere  in  space  and  vanishing 
at  infinity  together  with  its  first  derivatives,  can  be  represented  as  the  sum 
of  an  irrotational  vector  field,  AJ,  ,  and  a  solenoidal  field,  /U^.  •  The 

representation  is  unique  except  for  a  vectorial  constant.  Thus, 

ytf  =•  S10,  +  A0  £  U) 

where 

V  X  /U)t  -  0  (2) 

and 

V  '  Afi  ~  0  (3) 

The  application  of  this  important  theorem  to  the  study  of  wing-jet  flow 
fields  is  not  precisely  correct.  The  steady-flow  field  about  any  body  which 
sustains  a  resultant  force  due  to  fluid  pressure  will  include  a  vortical  wake, 
which,  along  with  the  velocity  field,  extends  indefinitely  if  the  fluid  is  inviscid. 
Such  flow  fields  apparently  do  not  meet  the  requirement  imposed  in  the  theorem, 
namely,  that  the  vector  vanish  at  infinity.  Nevertheless,  an  important  insight 
into  flow  problems  involving  free  vorticity,  as  is  characteristic  of  jet  flows, 
can  be  attained  from  the  following  considerations,-  based  on  the  more  restricted 
vector  fields  for  which  the  theorem  strictly  applies. 
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Lamellar  Vector  Fields 

It  follows  from  Eq,  (2)  that  AT/  can  be  represented  by  a  scalar  potential, 
such  that 


yZfj  =  —  +  const,  vector 

and  hence  satisfies  the  inhomogeneous  potential  equation 

V 2  (f)  ~  -  V  '  -  —  V  -  rur 


(4) 


(5) 


Consider  Green's  Second  Formula 


J(M  VzN-N  V2-M)c/r=J'('M^L-- N  (6) 

where  M  and  A/  are  two  scalar  point  functions,  continuous  in  the  first 
partial  derivatives  and  possessing  second  partial  derivatives,  and  /2  is  the 


outward  normal  seen  from  the  region  of  integration.  Let 


t?) 


N"~k 


(7) 


where 


yi  - 


TV 


M- fit' 


In  rectangular  coordinates 

Ji -'j(x-i)z+  +  (*-$)*■ 

where  %  t  /ij,  ,  -ir  are  the  coordinates  of  point  Tt  ,  and 

are  the  coordinates  of  point  M,1  .  So  defined,  M  is  harmonic,  i.  e.  ,  satisfies 

V 2  N  =  0 

everywhere  except  at  Jt  =0  and  hence,  except  at  this  point,  satisfies  all  the 
necessary  differentiability  conditions  everywhere.  The  volume  integral  in  Eq.  (6) 
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is  taken  in  the  Jtt'  system  over  the  region  enclosed  between  a  very  large 
sphere  of  radius  R  and  a  very  small  sphere  of  radius  €  >  with  centers  at 

the  point  JZ  =0,  Now  let  A//  be  identified  ar>  ^  of  Eq.  (5)  subject  to  all 
the  necessary  differentiability  conditions. 

Equation  (6)  then  becomes,  on  the  basis  of  Eqs.  (5)  and  (7), 


The  second  integral  on  the  right  of  Eq.  (8)  reduces  to  47T (j)  by  virtue  of  the 
assumed  continuity  of  0  .  Hence,  Eq.  (8)  can  be  written 

t TP  {-4r(v^r-  £  f  _  -i-  ff-]^  (9I 


Equation  (9)  is  an  expression  for  the  potential  at  points  within  R  in  terms  of 
the  divergence  distribution  and  the  values  of  both  and  its  normal  derivative 

on^the  boundary. 

* 

For  our  purpose  we  now  specify  the  vector  field  /Z/  to  be  such  that  the 
region,  or  all  the  regions,  in  which  V  •  /U  is  not  zero,  are  finite  in  extent 
and  can  be  enclosed  within  a  surface  Si 

In  Equation  (9)  consider  R  to  be  large  enough  so  that  if  R  j  is  the 
radius  of  the  smallest  sphere  that  encloses  Sf  ,  then 

R  »  Ri 

In  the  region  between  the  sphere s_^?_„  and 

V  ’  -  0 
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and  therefore  (p  in  this  region  satisfies 

V*  ^  =  o 

Since  (j)  represents  the  potential  throughout  space  of  a  suitably  restricted 
spacial  source  distribution,  it  will  vanish  at  least  as 

~R=R,~ 

This  situation  assures  that  the  second  integral  on  the  right-hand  side  of  Eq.  (9) 
vanishes  for  sufficiently  large  R  .  It  should  be  noted  that  the  corresponding 
velocity  field  also  vanishes  at  infinity.  Thus, 

<>«> 

It  is  shown  in  Ref.  17  that  ^  ,  as  given  by  Eq.  (10),  and  its  first  derivatives 

are  in  fact  both  defined  and  continuous  everywhere  in  the  field  if  •  /}/  is 
bounded  in  extent,  as  was  assumed,  and  is  at  least  piecewise  continuous.  The 
higher  derivatives  are  also  shown  to  exist  provided  the  divergence  distribution 
satisfies  a  Holder  condition. 


Solenoidal  Vector  Fields 


Let 


Mz  =  f  const,  vector 


with  the  condition 

V*  A  =  0 


(ii) 


(12) 
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By  forming  the  curl  of  Eq.  (11)  and  applying  the  vector  identity 

0 

VxVxJk=-V2JA  +  V(v-Jh)  (13) 

for  which  the  extreme  right  term  vanishes  by  Eq.  (12),  the  following  is 
obtained 

VZ  A=  “  =  VX/W  (14) 

Hence,  like  (fi  the  components  of  satisfy  the  inhomogeneous  potential 

equation. 

Again,  if  the  region  designated  ,  in  which  V  X  /yr  is  not  zero, 

is  finite  in  extent,  then  it  can  be  shown  on  the  basis  of  Green's  Second  Formula 
and  the  harmonicity  of  the  components  of  ^  exterior  to  5  £  that 

^  *  -jfp  J  fax  A s')  dt  (15) 

From  Eq.  (15)  one  may  observe  that  each  vector  element  of  rotation,  V X  SU?  , 
contributes  a  vector  element  parallel  to  .  However,  it  can  be  shown, 

since 

V-(V*/1*)=0  (16) 

that  the  net  effect  of  all  radial  components  of  V  X  ajj  (that  is,  radial  from  the 
position  of  A  )  vanish.  Hence,  more  concisely  it  is  the  tangential  component 
of  rotation-element  that  contributes  a  net  vector  element  parallel  to  A  . 

Evidently  all  the  conditions  concerning  the  necessary  functional  continuity 
of  the  distribution  V  ’  A/>  apply  to  the  components  of  \/x/ur  as  well,  and  in 
addition  the  condition  imposed  by  Eq,  (16)  must  also  be  met.  The  latter  condition 
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implies  that  the  rotation  distribution  itself  must  comprise  a  solenoidal  field. 

Further  consideration  will  be  given  this  fact  later. 

© 

Equation  (1)  can  be  expressed  in  terms  of  spacial  distributions  of  the 

0 

divergence  and  the  rotation  on  the  basis  of  Eqs.  (4),  (10),  (11),  and  (15). 

nj>  -  -  _ 1—  V clr+-  .1  v X  dr  +  vector  const  (17) 

4-TT  J  n  J-TT  J  Jb 

=  &+  ^tor 

where 

JTD.  =  V  X  AM 

For  any  given  vector  field,'  /l*  ,  which  satisfies  the  conditions  of  the  theorem, 

the  distributions  V'/Z 0  and  XlL  are  furthermore  unique.  (See  Ref.  18:  Chapter 
IV ;  Section  20.  ) 

Solenoidal  Flow  Fields 

In  steady-flow  problems  only  certain  information  concerning  the  flow 
field  is  given;  e.  g.  ,  the  conditions  far  upstream,  the  geometry  of  solid  boundaries, 
■  and  the  trailing-edge  conditions.  The  flow  field  must  be  determined  from  these 
data  alone,  and  its  establishment  involves,  in  essence,  the  determination  of 
appropriate  distributions  of  the  divergence  and  the  rotation.  Very  little  is 
known  concerning  distributions,  however,  except  on  surfaces.  Equation  (17) 
offers  a  means  by  which  surface  distributions  can  be  established  for  certain 
practical  problems. 

As  pertains  to  steady-flow  fields  of  an  inviscid,  incompressible  fluid, 

Eq.  (17)  requires  some  further  consideration  in  view  of  the  conditions  imposed  on 

Q 


© 
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the  vector  field  in  connection  with  the  derivation. 

First,  the  distribution  of  rotation  in  the  steady  inviscid  flow  field  about 
a  lifting  body,  in  contradistinction  to  the  condition  imposed  in  the  preceding 
development,  extends  indefinitely  in  the  wake.  For  our  purpose,  Eq.  (15)  is 
presumed  to  apply  generally,  provided,  of  course,  the  integral  exists. 

Second,  since  the  fluid  is  taken  to  be  incompressible,  the  entity  of  flow 
divergence  never  enters  into  the  actual  field  of  flow.  Rotation,  on  the  other 
hand,  will  certainly  be  present  in  the  wake  if  there  exists  a  resultant  force 
on  the  closed  stream  surface  of  a  body.  Provided  the  fluid  is  incompressible, 
the  former  statement  permits  us  to  formulate  any  steady-flow  field  which 
possesses  a  continuous  velocity  vector  (hence  one  which  vanishes  on  the  body 
surface)  on  the  basis  of  solenoidal  field  concepts  alone. 

On  the  other  hand,  if  the  fluid  is  also  inviscid,  in  which  case  the  velocity 
will  not  vanish  at  the  body  surface,  then  the  flow  field  must  be  continued 
analytically  into  the  region  within  the  body  surface  in  order  that  it  meet  the 
requirement  of  a  continuous  vector  field.  Under  such  conditions  it  may  be 
observed  that  the  flow  field  as  a  whole  does  not  comprise  a  definite  vector 
field.  But  this  situation  is  precisely  the  one  which  again  permits  us  to  treat 
such  flow  problems  on  the  basis  of  purely  solenoidal  flows  by  so  defining  the 
analytically  continued  flow  field. 

Since  we  will  be  dealing  exclusively  with  solenoidal  flows,  it  is  essential 
to  recall  that  the  rotation  distribution  is  subject  to  certain  basic  constraints. 
First,  the  divergence  of  the  rotation  is  everywhere  identically  zero 

v(vx^)  =  o 
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Therefore,  the  rotation  vector  field  itself  comprises  a  solenoidal  field.  The 
same  can  be  said  of  the  rotation  of  the  rotation  vector,  etc.  Second,  rotation 
is  a  convective  property  of  the  fluid  particles  subject  to  a  dynamical  requirement 
which  can  be  derived  from  Newton's  second  law.  Thus, 

JQJL.  - L —  Vj0+  (19) 

D  t  VP  ^  /o 

where  is  the  flow  velocity  (fyoo  +  /U>') 

-pj  is  pressure 

F  is  extraneous  force  (per  unit  volume) 

/O  is  density 

By  operating  with  the  curl  on  Eq.  (19)  and  applying  Eq.  (18),  the  following 

*  4> 

form  of  the  Helmholtz  equation  is  obtained 

where 

fl  =  Vxy 

*  • 

For  steady  conditions  the  preceding  equation  becomes 


Cv*v)/n=  (n-  v)  w+ 


(20) 


Equation  (20)  governs  the  rotation  distribution.  A  few  elementary  solutions  of  this 
nonlinear  system  have  been  determined  for  the  condition  of  force-free  flow: 
one  relates  to  a  circular  jet  in  a  parallel  stream.  Another  obvious  solution 
for  trailing  vorticity  is 

jflL  =  const  •  V 
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On  the  other  hand,  integration  of  this  equation  along  an  arbitrary  streamline 
leads  to  a  relationship  between  the  normal  and  tangential  component  of  the 

rotation  vector  which  will  prove  to  be  of  considerable  value  in  determining 

© 

the  vorticity  distribution  in  the  vortical  layer  of  a  deflected  jet. 

Consider  the  line  integral  of  Eq.  (20)  along  a  streamline  whose  vector 
element  is  designated 


J[(v-v)riJ-d*=f[(n-v)y].j*  +  / 


Since  \/  is  parallel  to 


as  follows 


f[(v-v)n]-djb 


A  ,  the  first  term  on  the  left  may  be  rearranged 

~fv-4^L  .  ctjb 
J  cLA 

~  f  V  Afii i  •  djb  +  /V 

*  CLA  J  CL  A 

=/  WdaT  -J  VC,  (A„  ■  C)  djL. 


where  JUj  and  are  respectively  the  tangential  and  normal  components  of 

M  referred  to  the  direction  of  the  streamline  element,  d.  .  The  operation 
involving  the  vector  derivatives  along  the  streamline,  which  lead  to  the  last 
expression,  is  carried  out  with  the  aid  of  Frenet's  formulas  (see  Ref.  19). 

Here,  C/  ,  is  the  First  Curvature  or  Flexure  of  the  streamline;  its  reciprocal 
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is  equal  to  the  first  radius  of  curvature,  ^  •  The  vector  C  is  a  unit  vector 

perpendicular  to  the  streamline  element  djb  and  in  the  direction  of  Jl  . 

The  first  term  on  the  right-hand  side  of  Eq.  (21)  can  be  developed  as 


follows : 


//(ti  =  /[(Sh-j-vyyJd.*! 

= fn.T l  v -pJfjriH ■  v)  vj.  cl*. 
=f/iTcCv+J'nN‘£L  . 
=Jflrdd+JnN-^  i* 


■  A> 


where  ctn^dni*  the  unit  vector  in  the  direction  Sh . 

The  preceding  developments  substituted  into  Eq.  (21)  give,  when  ff 
is  taken  to  be  zero 

fVclflT- fvC,SlN(in-c)clA=JflTdV+  J  Sl  N  cLV^  d. 


■  Ad 


or 


Jd&)  fvt  77  +  ("-V  d " 


Therefore 


Sr  + c,  v'  (*  ■  c)j + consd 


(22) 


Equation  (22)  is  an  expression  for  the  tangential  component  of  the  rotation  along 
a  streamline,  in  terms  of  the  integral  along  the  streamline  of  the  normal  component 
of  rotation.  The  tangential  rotation,  SI  f  »  includes  both  the  initial  value  of  this 
component  and  the  accumulation  of  the  "secondary"  rotation  which  is  known  to 


© 
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than  zero. 


III.  DEVELOPMENT  OF  THE  WING-JET  LIFTING  SURFACE  EQUATIONS 

r> 

We  now  apply  the  preceding  developments  to  the  formulation  of  the  wing- 
jet  problem  for  a  uniform,  circular  jet  under  the  conditions  of  small  deflections. 
The  integral  relationship  for  the  disturbance  flow  field  is  given  by  the 
solenoidal  part  of  Eq.  (17) 


where 

and  &LW  are  points  on  the  wing  surface.  Equation  (25)  can  be  arranged  so  that 
the  integral  is  taken  over  three  regions  where  the  rotation  is  known  to  exist;  the 
wing  region,  the  wing  wake  region,  and  the  jet  boundary  region.  Thus 
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cC^dycCQ 


l  K  dJJ  l [(*-&¥+  -£?]*■„ 


JjJ t[(*-W+(r^+(^Y]yt jd*' v 


(30) 


fS\ 


If  the  local  surface  inclination,  is  everywhere  small,  then  the  wing 

and  wing  wake  terms  can  be  simplified  and  combined  to  give  the  usual  surface 
integral  relationship  for  downwash  from  lifting  surface  theory.  Thus,  for  the 
wing  and  wing  wake  terms,  the  first  triple  integral  term  on  the  right-hand  side 
of  Eq.  (30)  can  be  neglected  in  comparison  with  the  second,  the  individaul  terms 


being  of  higher  order 


in  f‘(y) 


wa/ 


than  those  of  the  second  triple  integral. 
Therefore,  the  first  two  terms  on  the  right-hand  side  of  Eq.  (26)  can 
be  represented  approximately  as  follows 


Vr-t-  vs- 

xv  J w/  ^  vv 


where 
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The  terms  on  the  right-hand  side  of  Eq.  (31)  can  be  rearranged  into  the  following 


form 


.k  J-  C 

i  2 


[&-%)•+  (^r]%  A  i.£ 


Partial  differentiation  of  the  first  term  above  with  respect  to  and  the  second 

term  with  respect  to  <£>  ,  together  with  the  divergence  relationship 


X  _ 


yields  the  following  expression  for  Eq.  (31) 


vs  j*+vs. 


k  J-  £■ 

2  2. 


-  /  \  u  (33) 


-b  [  //V (~p'oi'6')- *V  ('r6 - *d] +  IXy- (*• + e)~ (*•’ e)X 'Ij.t 

J-jl  [  (*-%)/  (*-$,)*- f  (y  +  4.)2-  C *-$) /&-&)*+ (f-fc) v 

It  should  be  stated  concerning  Eq.  (31)  that  we  anticipate  a  discontinuity 
in  Ajs  af  the  wing-jet  boundary.  Equation  (32)  therefore  does  not  apply  at 
such  points  and  Eq.  (33)  is  developed  in  that  light. 

Consider  now  the  third  term  on  the  right-hand  side  of  Eq.  (26).  Its  form 
as  given  by  Eq.  (30)  can  be  expressed  in  cylindrical  coordinates  as  follows, 
with  the  x-axis  of  the  two  systems  coinciding. 
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y  =  n  cos  z^ 
-z-  =  /t  s/n  i? 


-  jn'  cas  y 

-  ft1  s/n  y> 


and  therefore 


i7z  =  fl-i/  cosy  +  fl^smy 
-ft  y,  =  -fi-V'Siny  +17^  cos  y 


vs  •  J^J-r  dr  =  f  fo  fff [(n 

Iff  L~(?  "4}  f^-iy^in^f-Sl/L.icosTp)  -  n(co$'lS'-n'  cts^Q*]  ,  ,  ,  i35^  ^ 

JJJ~ldtr~* + o&^^s  *-'*•'?*  * 


Equation  (35)  will  be  simplified  under  the  conditions  that  the  jet  is 
initially  circular  and  undergoes  only  small  deflection  and  cross-sectional 
distortion.  The  jet  rotation,  Jfk  j  ,  will  in  that  case  be  essentially  in  the  region 
A/  equals  a  constant,  Rc  ,  and  terms  containing  -fl-n.  can  be  neglected. 


We  may  write 


~S'%y'  -  J'SLy  cL  ft1 

Yx  -  J fix  dnJ 


On  this  basis,  the  expression,  for  the  downwash  on  the  wing  surface,  due  to  the 
jet  vorticity,  becomes 
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C  zir 

vs-  ir-.Rcf(x)  /  /  - 

J  «iJ  J„c/0  ffr-Qy+f*-!-  JF-  %  «s  <fl  % 

+%  r  f **&' ^  ^ s'«y) -(<t- r„  as y)  rfj  dxfjt  (3 

J-coJo  'll2*  RS-2y  R„  cslfjvi  ™  * 


where  the  wing  surface  is  taken  to  lie  in  the  plane  'if' =0. 

To  a  first  approximation,  ,  in  the  first  term  on  the  right-hand 

side  of  Eq.  (37)  may  be  taken  to  be  the  constant  value  of  jet  vorticity  far  upstream, 
Y  ,  in  which  case  the  expression  can  be  integrated  at  once  to  give 

Wo  s  \ 


Astor-t,  - _ 


,,  y  /  '47r  -t"  <  x,~) 

-  -/Ml,  +  ilf  M  y  ^  R.  I  , 

•  (_  0  £°'r  ->  Re  J 

The  other  surface  integral  on  the  right-hand  side  of  Eq.  (37)  can  be 
developed  by  integration  by  parts  allowing  for  the  fact  that  Yqy  will  be  discon¬ 
tinuous  at  the  wing-jet  boundary.  Thus,  partial  integration  of  the  T^xterm 
with  respect  to  ~lj)  yields 


[+°°  f  $>)Yv*sm  y>  cttfcSS, 

Jo  r  (Y-tY-h^  +  R*-2^R„ 


Jo  [  i-  +  R0Z-  RecosZf]  V2- 

=  _[$(><'$,)  ffr(vrt)’7Uw-0  , 

J-  s.  ^  £  \/ (x-  £>)*'  ■f'fyt’R*)*  | /7Y-^yTT*pR<> 

•  e  P“  (x.-&)  f _ £ _ 


-y  d£> 


^  £  fa- **‘7 
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where  use  has  been  mad^  of  the  fact  that 

ty.  (rr-e)  -  r^(t) 
7, *  (ift-  e)  =  r(rr+i) 


/  jYl 


'ly 


and  also  of  the  condition  4  j 

r*x  _ 

>6  /?.  ■3? 

at  appropriate  points. 

Integration  by  parts  of  the  Yx  term  on  the  variable  £? 

’ +oo  f>2,TT> 


4 


-  (y-Ko  cestp)Yx  dy>c/ & 


/-  tfO 


= 


/V"  +*S-Lr*o*sv] 


?Yx 


r+°° ri7r  (y-Ka  <*sv)(x-e,)-yi-cLzfJG> 


gives 


*>S  ?P)  Y/C°°)  c/ if 


(40) 


'~a> 


The  total  contributions  of  the  jet  vorticity  to  downwash  at  the  wing  surface, 
as  given  by  Eq.  (37),  consists  of  the  sum  of  Eqs.  (38),  (39),  (40).  These  relation¬ 
ships,  together  with  the  wing  vorticity  effect  as  given  by  Eq.  (33)  and  also 
Eq.  (29),  substituted  into  Eq.  (26)  yield  the  following  expression  for  the  total 
downwash  on  the  lifting  surface. 


■v, 
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^  J 

ri 

1  n  1 

>  K 4 

2-/U-  ca  s 


y*+x0t-*Y*.*'*y J 


j_  j  _f  Alfa  (~Rq)  ('ft+R*) _ -/-  (R«)  ('Y  ~  _ 

H  l  ft- tty#-#*’ +(?+**)*  (*'§>))/  (x-t,) *  i{$'Rtr 

-  a'C'M  (*-£)  _  -Af^fc)  (x-t*)  ) 

^ /<" +  fy+x*)'2'  i / ’*- z ^ ^ ^ 4  J 

where  /  \  , 

=  '!?(-**+£)  -  7y(-R6-e) 

AK'iy  (+R<7)  —  /y  (-t  Ratf)  —  Yy,  (+£„-€.) 

( rt)  -  7^  (rr+i:)  ~~  7^  (rr- {) 

^fa(o)  ~  1^(0  ~  Yvsfzrr-c) 

The  jet  trailing  vorticity  distribution  designated  Y %  in  the  above 

equations  can  be  determined  approximately  by  the  use  of  Eq.  (22),  which  applies 
along  streamlines  within  the  jet  vortical  layer.  Along  any  particular  streamline 
not  subject  to  body  forces,  we  have 


+  (R„)  fy-/?J 

'*  (X'$?y (x-iY  tty'*,) 

-  AiirCo)  (x-t*)  ) 

^  /( X'  fy) 1  +-  fa  -  ^  *  f 


where 


YY  r  — 
V 


=r^[-£Y+<'v(”-c)j^ 
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It  can  be  shown  that,  in  the  absence  of  body  forces,  the  velocity  derivative 
in  the  direction  of  ft  ,  which  is  both  normal  to  the  streamline  and  tangent 
to  the  Bernoulli  surface  is  also  equal  to  V •  Ctt  where  d,  (m  •  <c  )  is 

the  geodesic  curvature  of  the  streamline  relative  to  its  Bernoulli  surface.  In 
Cartesian  coordinates,  the  flexure,  Q f  ,  may  be  expressed  as  follows 


For  streamlines  curved  only  slightly  from  their  original  X  direction, 
d.A’  is  approximately  cLX  ,  and  the  above  expression  can  be  represented  by 


4  = 


dir'2' 


'd2-* 


Furthermore,  for  small  values  of  jet  deflection  and  cross-sectional  distortion, 
(tl  •  C  can  be  evaluated  on  the  cylinder  ft0  -cots/- in  which  case  it  can  be 


expressed  as  follows 


Pi  '  €  -  CCS jfjZ-  -  TArf 


/  cl** 

'/  f  d 

-4-^ 


where  ^  is  the  azimuthal  coordinate  introduced  previously. 

Under  these  conditions,  Eq.  (42)  may  be  written  as  follows 

=  2  f  ^  f-stn  Zp  jillt  +.  C09  ZO  dl±-  7dx  (45) 

v  J.oo  v L  1  dxl  >  dz^J 

If,  in  addition,  V  and  dl ^  are  assumed  not  to  vary  mufch  along  streamlines, 
the  integration  can  be  performed  and  one  obtains 
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nT-  ZflH  / -s/n  If  ^  -f-  cos  ZJ>  jL 
'  J  diC  '  cL 


cL  ^ 


(46) 


Integration  of  Eq.  (46)  across  the  jet  vortical  layer  yields  an  expression 
for  the  secondary  vorticity.  Here,  we  apply  the  Second  Theorem  of  the  mean 
for  integrals  to  obtain 


v- 2  v  /- sin  ?#• + «"  y-gf 


cL, 


(47) 


where  and  are  appropriate  mean  values  of  the  strean*line 

deflections  across  the  vortical  layer. 

The  Y*  component  of  the  jet  vorticity  can  now  be  determined  for  the 
condition  of  small  deflections  by  taking  appropriate  components  of  %~Y t  and 
!  the  result  is 


h  -  tfjz-  *  -  vfr] 

It  may  be  noted  that  the  magnitude  of  the  jet  vortical  layer 


(48) 


depends  on  the  product  of 


X*.  =  v,-v„ 


J  Yon 


and  the  mean  deflections  of  the  stream¬ 


lines  which  are  of  order  .  For  small  values  of  \J j  -  ,  i.  e.  of 

order  /  *(x)  ,  Yt"  ts  therefore  of  higher  order  than  the  wing 

vorticity  terms  retained  in  Eq.  (41)  and  accordingly  such  terms  can  be  neglected 
in  determining  first-order  effects  on  the  wing.  On  the  other  hand,  the  lift 
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distribution,  to  a  first  order  in  small  values  of  hZ. — ,  as  for  example  Koning, 

Vx> 

Ref.  (1),  must  necessarily  be  second  order  in  terms  of  wing  angle  -  of-attack. 

© 


© 


The  mean  deflections  and  dz/fCX  may  be  expressed  in 

terms  of  the  normal  disturbance  velocities  evaluated  on  the  jet  vortex  sheet. 
In  that  caseT  *  * 


^  _  _ 

Voo-hVj  ^  Voo  +■  Vf 


Tfer  _  2.  JjJ  -h~UJ~T 

•  LU  J 


(49) 


where  /i^  and  are  due  to  the  jet  vorticity,  and  and  are  due 

to  the  wing  vortex^  system. 


**  =- 


ztr  =  —  -L 


1  OO 


L  £  ^  -  C 

2  2 


Substitution  of  Eqs.  (49)  and  the  first  two  of  Eqs.  (50)  into  Eq.  (48)  yields 
the  following  expression 


?.4 


e 
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I 

. 


i0°  raff  )]YX(&,  Zf)drfd& 


rr*LL _ _ _ ___ 

m(Voo+Vr)  Jo  £(x-t,y+  -  C0s(siy-V})jj3/*. 

%-A)  stn  C^-y)  d^pdk, 

[(X '$)*+**.*[  -  cos  4'i- 


-r, 


V^f_ 


-too  rzir®/ 


(51) 


'L  ^  'Kj  sm  ir  +  2-Y^  KS  00 s  7^ 

_ * _ _ _  _ ^ _ _ _ 

/  /  ^ 


V#>  -h  Zj 


Yoo  +■  Zt 


Let 


K  =  Y*  (xtf) 

x?ty 

V  (Si1?)  -  r^*^  +  l*j£  (V*)  +  ~"' 

Since  Vy  (X,£?)is  symmetrical  about  the  jet  periphery,  the  contribution  of  the  first 
term  on  the  right  hand  side  of  Eq.  (51)  depends  on  the  higher  differential  coefficients 
of  the  Tavlor  expansion  than  the  first,  and  therefore,  according  to  Eq.  (48)  on 
the  derivatives  of  the  streamline  deflection  higher  than  curvature.  The  contribution 
of  the  second  term,  however,  depends  on  the  first  differential  coefficient  of  the 

expans  ion.  It  may  be  noted  that  even  this  coefficient  already  depends  on 
the  product  of  streamline  deflection  and  curvature.  For  small  values  of  deflection, 
curvature,  etc.  ,  Eq.  (51)  may  therefore  be  approximated  by 


V,  = 


tvZj  sm  2^ 
.  Vx  +■  Zoo 


+ 


cas'lt' 

Zoo  ZT 


(52) 
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The  last  two  of  Eqs.  (50)  can  be  expressed  as  follows 


O 


2 


<ur 


^  _  R*  vniy _ 


/  ~  37^  / 

T 


/  + 


6r-Q  )./iA 

/?2r-^;v  2rW  d*>dyi 


JllT" 


- _ / 


c 

'L 


/  i  fc-fi) _ _  //%-Kj  0*-re'.7j) 

irr  J_L  [(x-k,t+R? «»*&/  £  focHy^ T^nT^y) 


^/~ 

[77TPT7[J7R- 
_ f_  C*  f*  avy  Ro  **  ly-n  f .  , 

*2  l 

(x^)[(x.kr  +  r+^s,^j  ? 

[(%-t,  Y  +  Zc’-S"’1  Vj/Oc^b1) 

where  b*' =  R0X  ~  Z  R0Y[  co s  V' +  rj'1 


(53) 


Upon  substitution  of  the  above  expressions  into  Eq.  (52)  one  obtains 


ATT 


1Y« 


*¥■ 


.  [Z  2P  (R0-n  cos'll)  / 

L  9v  S2  ? 

a 

r [(*%&*]* 


(54) 


26 


AI-1190-A-6 


e 


t 


.  (l3S*o  COS  77  cos'll) 

'  (*-  +  Roz  *,n%  w 


cLfey  cLy^ 


Ro_  cos 

4-rr 


,  ©  L  \ 

(Y-i) 


fry,  (-Rt+€ )- >V  (-R.-C  )/  0+ccj  tT) 


/- c. 

2 


(V-W  +  RSS-'V*  £  ]/(r-t9ytiRS(ti-cos*>') 


__  fry  (*•+£)  't*  (*t -£)](/'L°sfr) 

y(x-i¥  +  2R*  (i-coslp) 

For  large  values  of  ^  ,  Eq.  (54)  becomes 

( V'-v-Zr  \  y,  /i.s'i  - - L.  [*  JT_  (k.:%™*L 

tTrH  5? 

and  the  derivative  of  Eq.  (54)  with  respect  to  %  yields 


•r 


(55) 


y  00  +  Vr  1  _ l_ 

1*  rr 


b 

1 


t 


(56) 


('Z’VJ  5*  n  J.tJ.c  3? 

'  2.  z. 

R  -Y)cosiy*+-  ^  (** 

0  ^  (Y'l,)2  +  R0l  St nl  if 

_  Zb2(t-b)L  (R„  cosz  'iT-r]  cos  fr) 

i fhW+V  [(x-bY+R.1  *•"'*] 

By  considering  once  more  the  rotational  flow  aspects  of  the  vortex 
sheets,  it  can  be  shown  with  Eq.  (20)  that  if  the  wing  lift  is  considered  to  be  a 
constant  across  the  jet  rotational  layer,  there  will  be  no  shed  vorticity  arising 
at  the  wing  and  jet-boundary  intersections.  Therefore 
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(57) 


A  fy  (+  ~~  A  Yv  Co) 

and  Ro  )  =  A  Tfy. 


Equations  (41),  (55),  (56),  and  (57)  complete  the  formulation  of  the 
lifting  surface  problem  for  a  wing-jet  combination  for  the  condition  of  small 
flow  deflection.  -Since  the  system  of  equations  appears  quite  formidable  in  its 
present  form,  one  must,  as  in  the  case  of  ordinary  lifting  surface  theory, 
resort  to  approximate  methods,  employing  expansions  that  are  suitable  for 
special  classes  of  configurations  in  order  to  obtain  numerical  results.  It 
appears  that  wing  aspect  ratio  offers  a  suitable  expansion  parameter  for  this 
purpose,  since  the  limiting  cases  of  both  high  and  low  aspect  ratio  yield  great 
simplification  in  ordinary  lifting  surface  theory  which  is  implicit  in  the  present 
p  roblem. 

This  course  has  been  pursued  for  the  case  of  a  high  aspect  ratio  wing  and 
is  presented  in  the  Appendix.  For  simplicity,  jet  curvature  terms  have  been 
neglected.  The  lifting  surface  equations  given  here  in  that  approximation, 
reduce  to  a  couplet  of  integral  equations,  the  same  as  given  by  Graham  et  al. 
in  Reference  13  but  with  somewhat  different  expressions  for  the  image  strength 
factors.  In  phce  of  the  quadratic  expressions  (i/j  -  VjJ/  (Vj-V  Vgo)  and 
ZVj-Vg, jC^^oa )  which  occur  in  basic  image  analysis  of  a  perturbed  circular  jet, 
the  present  development  obtains  linear  expressions  (Vj-~  Vj  and 

ZVo/VrtVco  respectively.  Numerical  results,  showing  the  effect  of  these 
differences  are  given  in  Figure  2  where  a  comparison  of  the  wing  lift  increment 
due  to  the  jet  is  made  with  other  existing  theories  and  experiment  for  one  set  of 
conditions.  The  present  developments  yield  values  of  wing  lift  within  the  jet  region 
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A  formulation  of  the  lifting  surface  equations  for  a  wing -jet  configuration 
is  given  for  the  case  of  a  circular  jet  experiencing  small  deflections.  Jet 
curvature  effects  are  included  but  jet  cross-sectional  distortion  is  neglected. 

It  is  shown  that  any  steady  flow  of  an  incompressible,  inviscid  fluid 
can  be  represented  as  a  solenoidal  field  with  the  corresponding  rotation  distribu¬ 
tion  subject  to  the  dynamical  equations  of  rotational  flow.  A  useful  streamline 
integral  of  the  latter  system  is  derived  that  gives  the  distribution  of  the  stream- 
wise  component  of  the  rotation  along  curved  streamlines  within  the  jet  rotational 
layer.  This  is  identical  with  Hawthorne's  expression  for  secondary  vorticity 
in  nonuniform  flows.  The  flow  field  of  the  wing-jet  combination  is  first  rep¬ 
resented  as  a  general  solenoidal  field.  It  is  then  specialized,  for  the  case  of 
a  lifting  surface  and  an  initially  uniform  circular  jet,  to  one  involving  only 
vortex  sheets. 

On  the  jet  sheet  both  the  normal  and  the  streamwise  components  of 
vorticity  enter  into  consideration.  Hovfever,  the  downwash  effects  of  the  jet 
on  the  lifting  surface,  with  exception  of  discontinuity  effects  at  the  wing  and  jet¬ 
boundary  intersections,  can  be  expressed  entirely  in  terms  of  the  trailing 
component  of  jet  vorticity  and  its  streamwise  derivative. 

With  the  aid  of  the  aforementioned  streamline  integral  of  the  rotational 
flow  equations,  an  approximate  expression  for  the  jet  trailing  vorticity  is 
obtained  in  terms  of  the  wing  vortex  system.  For  small  values,  the  jet  trailing 
vorticity  depends  on  the  product  of  the  initial  jet  vorticity  far  upstream  and  that 
component  of  the  local  flow  deflection  which  occurs  due  to  geodesic  curvature 
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of  the  streamlines  on  the  jet  vortex  sheet.  Its  value  is  a  maximum  on  the  sides 
of  a  downward  deflected  jet  and  vanishes  at  top  and  bottom  with  sign  such  as 
to  cause  upwash  within  the  jet  and  downwash  on  the  outside. 

For  simplicity,  jet  curvature  effects  have  been  neglected  in  that  develop¬ 
ment.  The  resulting  system  then  reduces  to  a  couplet  of  integral  equations  the 
same  as  given  by  Graham  etc.  in  Reference  13  except  with  image  strength  factors 
altered.  Numerical  results,  showing  the  effect  of  these  differences,  yield 
values  of  wing  lift  within  the  jet  region  considerably  higher  than  previous 
theories  and  available  test  data  for  a  nonuniform  jet. 

The  approach  taken  here',  and  basic  developments,  lend  themselves  to  the 
investigation  of  other  problems  involving  jet  curvature  effects.  Thus  whenever 
a  jet  is  turned  toward  alignment  by  a  main  flow,  there  is  as  a  consequence, 
and  in  accordance  with  the  general  prescription,  a  redistribution  of  its  vorticity, 
the  immediate  effect  of  which  is  a  self-consistent  flow  field  both  within  and 
exterior  to  the  jet.  Through  the  extended  influence  of  this  flow  field,  jet 
curvature  effects  manifest  themselves  on  solid  boundaries  of  wings,  propellers, 
and  other  similar  parts. 
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APPENDIX 


*»  LIFTING  LINE  APPROXIMATION 


.  If  the  high  aspect  ratio  approximation,  i.  e.  ,  |  j  |  ^  |  ’ 

which  ordinary  lifting  surface  theory  is  reduced  to  lifting  line  theory,  is  adopted 
and  applied  to  Eq.  (41)  with  the  jet  curvature  effects  term  neglected,  an  integral 
equation  is  obtained  for  the  spanwise  circulation  distribution  from  which  numerical 
results  can  be  derived  quite  readily?  Under  such  conditions,  Eq.  (41)  with 
given  by  Eq.  (54),  becomes 


trfjOlV"  +) 


Ro_  Vj  -  Vm 

*  Yr^*> 


vro,  *+<  R0 

0  f  ^  Rg 
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The  above  equation,  for  the  case  of  a  flat  plate  wing,  i.  e.  ,  (x)  =  oC  ,  can 


be  written  as  follows 


where 


(2-A) 
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Inversion  of  Eq.  (2-A)  yields 
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which,  with  as  given  by  Eq.  (3-A),  becomes,  upon  integration 

r(ft)--zg.  Ffa) 


(4-A) 


Equation  (4-A)  with  as  defined  following  Eq.  (2-A),  is  an  integral  equation 

for  the  wing  spanwise  circulation  distribution,  /^(^)  . 

In  the  expression  for  F (ty)  the  term  which  is  a  surface  integral  can  be 
reduced  to  a  single  integral  by  performing  the  integration  on  .  Equation 

(4-A)  becomes  therefore  for  £?<f  ^  <£  Rc  * 
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Equations  (5-A)  and  (6-A)  are  equivalent  to  Eqs.  (3,  11)  of  Reference  13,  which 
were  developed  on  the  basis  of  simple  image  theory.  In  fact,  if  the  present  system 
is  integrated  partially  one  obtains  Eqs.  (3,  11),  but  with  different  expressions  for 
the  image  strength  factors  that  precede  the  integrals.  Instead  of  ~ 
and  2  VjV^  j (y j  +^0  ^  the  present  expressions  turn  out  to  be  (V] j  (y 

and  W„/(Vj+V„) 

The  procedure  used  here  for  solving  Eqs.  (5-A)  and  (6-A)  makes  use  of 
a  single  sine  series  expansion  far  I-1.  The  leading  terms  of  the  series  were 
determined  in  the  usual  manner  by  solving  sm,  simultaneous  equations  which  apply 
at  ori  discreet  points  along  the  semi- span  of  the  lifting  line.  For  this  purpose  it 
was  convenient  to  let  H  be  expressed  as  follows 


r,('f)=TTCcLVTLAt^i  Q 


The  section  lifts  coefficient,  "defined  as 

c  -  y'v r 

‘  *  J£S0VJC 
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can  then  be  written  as  follows 


(i„.  ,), 


where 


)Vj"  inside  the  jet 
\4o  outside  the  jet 


Numerical  results  for  the  spanwise  lift  distribution  obtained  with  a  five 
point  collocation  solution  of  the  present  systems  of  equations  are  shown  in  Fig.  2. 
For  these  computations,  the  five  points  were  taken  according  to  the  following 
formula 


=  s  cos  <n- 


/n  =  1,2,  3,  4,  5 


and  the  remaining  parameters  were 


=  .  735 


Wing  aspect  ratio  =  5 

Z$°  =  .  12 
b 


The  theoretical  curves  in  Fig.  2  show,  in  essence,  a  comparison  of  the  lift 
increments  due  to  the  jet  as  computed  by  the  various  existing  theories.  In  the 
present  analysis,  this  increment  was  computed  by  subtracting  out  the  theoretical 
lift  for  the  wing  in  a  uniform  flow*obtained  by  using  Glauert's  lifting  line  theory. 

The  unusually  high  value  for  the  maximum  lift  predicted  with  the  present 
theory  as  compared  with  that  of  Reference  13  reflects  the  effect  of  the  difference 

in  the  constants  cited  above,  as  well  as  a  possible  difference  which  may  arise 

£> 

from  the  respective  methods  of  treating  the  equation.  In  Reference  13  two  sine 
series  were  employed  in  the  expansion  of  r  ,  one  applying  outside  the  jet  and 
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the  other  inside  the  jet  with  a  matching  of  the  lift  at  the  wing®jet  boundary.  A 

single  sine  series  expansion  was  assumed  for  /"*  in  the  present  computations. 

It  should  be  noted  further  that  the  present  lifting  line  analysis  neglects 

jet  curvature  effects,  and,  in  addition,  is  based  on  the  condition  that  Vj-  is 

0  / 

not  close  to  ,  more  specifically,  that  V qq  is  of  order  Veo 

The  numerical  computations,  however,  have  been  carried  out  for 


Voo  =-36>0 

which  is  close  to  the  range  where  wing  effects  to  the  second  order  in  angle  of  attack 
may  be  of  importance. 

In  view  of  seemingly  poor  agreement  with  experimental  data,  it  is  at  this 
time  difficult  to  justify  the  present  method  over  previous  analyses  of  the  wing- 
jet  interaction  problems.  Nevertheless,  it  is  felt  that  the  present  approach, 
with  the  jet  boundary,  regarded  as  a  free  vortex  sheet,  affords  a  more  realistic 
model  than  does  the  image  technique,  as  well  as  a  better  understanding  of  the 
effects  of  the  controlling  parameters,  i.  e.  ,  jet  deflection,  curvature,  etc. 

If  the  results  of  simple  image  theory  of  Reference  13  are  to  be  recovered, 

t 

one  is  led  to  conclude  that  either  Eq.  (49),  or  the  ensuing  developments  leading 

to  Eq.  (52),  must  be  modified.  This  follows  from  the  fact  that  the  expression 

for  on  the  jet  boundary  as  given  by  Eq.  (48)  corresponds  exactly  to  that 

which  is  obtained  with  the  flow  fields  from  simple  image  analysis  if  and 

.fiL  are  taken  to  be  the  mean  deflections  of  the  streamlines  on  the  two  sides 
(Lnii 

of  an  elementary  jet  boundary.  Equation  (52)  on  the  other  hand  cannot  be  so 
reconciled.  The  need  for  further  analysis  is  suggested  in  this  area.  In  addition, 
it  is  recommended  that  a  more  complete  development  in  the  lifting  line  approximation 
be  carried  out  to  include  the  jet  curvature  effectf  terms. 
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It  should  be  stated  that  experimental  data,  more  consistent  with  the 
theoretical  flow  model  than  is  at  present  available,  is  required  in  the  final  correlation. 

•  © 

The  ideal  jet,  which  is  usually  assumed  for  simplicity  of  analysis,  can  perhaps  best 

f?) 

be  simulated  in  a  wind  tunnel  by  employing  drag  screens  rather  than  a  propeller 
whose  wake,  due  to  its  characteristic  loading,  is  generally  nonuniform.  This 
is,  of  course,  not  to  lose  track  of  the  ultimate  problem  which  pertains  to  a 
real  propeller  and  wing  combination,  but  merely  to  test  the  ideal  jet  model.  It 
could  well  be  that  the  realistic  problem  requires  an  even  more  elaborate  model  for 
the  propeller  wake  such  as,  for  example,  a  s*eries  of  concentric  ideal  jets. 
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WING-SLIPSTREAM  ARRANGEMENT  SHOWING  COORDINATE  SYSTEMS 
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